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ABSTRACT 

Using the outer gap model, we investigate the emission region for the multi- 
wavelength light curve from energetic pulsars. We assume that 7-ray and non-thermal 
X-ray photons are emitted from a particle acceleration region in the outer magneto- 
sphere, and UV/optical photons originate above that region. We assume that 7-rays are 
radiated only by outwardly moving particles, whereas the other photons are produced 
by particles moving inward and outward. We parameterize the altitude of the emission 
region as the deviation from the rotating dipole in vacuum and determine it from the 
observed multi-wavelength pulse profile using the observationally constrained magnetic 
dipole inclination angle and viewing angle of the pulsars. We find that the outer gap 
model can explain the multi- wavelength pulse behavior by a simple distribution of emis- 
sivity, and discuss the possibility of further improvement. From observational fitting, 
we also find a general tendency for the altitude of the 7-ray emission region to depend 
on the inclination angle. In particular, the emission region for low inclination angle is 
required to be located in very low altitude, which corresponds to the inner region within 
the last-open field line of rotating dipole in vacuum. This model suggests a modification 
of statistics about observed 7-ray pulsars. Number of the sources with low inclination 
and viewing angles increases compared with previous estimate. 

Subject headings: stars: magnetic field — stars: neutron — pulsars: general 

1. INTRODUCTION 

Pulsars emit over a wide range of energies from radio to 7-ray. Recent observations by the Fermi 



Gam ma-Ray Space Telescope of more than sixty pulsars (jAbdo et al.ll2010al : ISaz Parkinson et al 



2010) have revealed further details of the structure of the emission region. The detection of the emis- 
sions in the GeV energy range from a pulsar magnetosphere means that electrons and positrons are 
accelerated to more than ~ 10 12 eV by the electric field paral lel to the magnetic field, which arises 



in a depleted region of the Goldreich- Julian charge density (jGoldreich Julian! 11969 ). The light 



curve in the 7-ray band is an important tool for probing the particle acceleration and dissipation 
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processes in the pulsar magneto-sphere, since the maximum energy is determined by the acceleration- 
radiation-reaction- limit for typical energetic pulsars. The 7-ray emission reg ion has therefo r e bee n 
explored by comparing theore t ical m o dels with the observed ligh t curvefe.g.. IWatters et al,l (|2009l ); 
Venter. Harding fc Guillemot! (|2009l ); iRomani fc Watterd (|2010l )). The pulsed emission is also de- 
tected in other energy bands (X-ray, ultraviolet, optical and radio) for some sources (e.g.. lThompson 
(|2004l )). The spectral features are non-thermal except for the soft X-ray range, and the light curves 
from a single object are, in general, different from one energy band to another. For example, profiles 
of the light curve i n one spin period are different in the 7- and X-ray ranges in the Vela pulsar 
(jAbdo et al.l l2009al ) . The peak phase of different energy range is expected to coincide, since the 
emitting particles are related to a pair cascade process. However, the observation shows that the 
phase depends on the energy bands. This means that their emission regions are not the same region. 
A complete understanding of light curve behavior in multi-wavelength bands can provide valuable 
information about the particle acceleration region. Note that we do not discuss soft X- r ays, w hich 
are believed to be thermal radiation from the neutron star surface (e.g., iJackson et al.l (|2002l )). 



Possible ori gins of non-thermal pulsed em issions have bee n considered in the polar cap (IDaugherty &; Harding 

1996 ). slot gap ( Muslimov Harding 2004 ). and outer gap ( Cheng. Ho Sz Ruderman 19861 ) mod- 
els. Recent Fermi observations with high 7-ray photon number statistics have showed that the 
phase-averaged spectrum above 200 MeV is well fitted by a power law plus exponential cutoff, 
and that a c utoff s hape sharper than a simple exponential is rejected with high significance (e.g., 
Abd o et al.l (12009ah). This ru les out the near-surface emission proposed in polar cap cascade models 
(jDaugherty fc Harding) 119961 ) . which would exhibit a much sharp spectral cutoff due to magnetic 
pair-production attenuation. Thus, pulsed 7-ray emission originates in the outer magnetosphere, 
as considered in the outer gap model. 



Takata. Chang Shibatal (|2008l ) (hereafter TCS08) considered a three-dimensional geometri- 



cal emission model to fit the observed light curves at different energy bands. The model is extended 
with some model parameters from the calculated gap structure in a two-dimensional meridian plane. 
By comparing the light curves of the Vela pulsar, they found that the X-ray emission is produced 
by both inward and outward emission from the gap region, and that UV/optical emission origi- 
nates from secondary pairs at a higher altitude. The number of light curves of pulsars observed at 
7-ray and other energy bands is increasing thanks to Fermi, so that it is worthwhile to investigate 
whether outer gap model is applicable to other sources. 

In this paper, we investigate the emission regions of several pulsars by fitting the simplified 
model of TCS08 to the observed multi- wavelength light curves. In this model, we have to specify the 
locations of the upper and lower boundaries of the gap region where the non-corotation potential 
is zero. Therefore, we explicitly introduce the altitude of the gap region as a parameter, in order to 
fit the observational data easily. The light curves also depend on the dipole inclination and viewing 
angles. In our method, such parameters are eliminated by other observational data, and only the 
altitude is changed for the fitting. In the most studies, the lower boundar y of the emission region is 
chosen as the surface of the last-open field lines of the rotating dipole(e.g. iTakata. Chang &; Shibata 
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(|2008l ); iRomani fc Watterd (|2010l )) In this paper, however, the altitude is allowed to be in a wide 
range in order to explore the possible deviation of magnetic field-line structure from that of a 
rotating dipole in vacuum. In Section 2, we describe the model assumptions and parameters. In 
Section 3, we compare the peaks of light curves with those observed at multiple wavelengths, and 
determine the altitude parameter. Our discussion is presented in Section 4, and lastly, a summary 
is given in Section 5. 



MODEL DESCRIPTION 



The numerical me t hod fo r fitting the light curve is well described bv lRomani fc Watterd (|2010l ) 
and lBai &: Spitkovskvl (|2010al ); but we briefly summarize it in this section to explain one modifica- 
tion. Our model is almost the same as that used by TCS08. However, we explicitly introduce the 
altitude of the emission region as an additional parameter. 

We assume that magnetic field structure is approximately described by a rotating dipole with 
magnetic moment /i. The angular velocity is O, and the magnetic axis is declined by an angle a 
from the axis of rotation (the z-axis). The magnetic moment changes with time t as 

= //(sin a cos ilix + sin a sin £lty + cos az). (1) 

The magnetic field produced by the rotating dipole (e.g. Jackson (jl975l )) can be expressed using 
the retarded time t r = t — r/c as 
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(2) 



where r is radial distance from the center of the star, and a dot denotes a derivative with respect 
to t. 

We assume that radiation direction aligns with magnetic field in a frame rotating with angular 
velocity fi in which the electric field vanishes. Physically, this means that the magnetosphere is 
filled with a co-rotation enforcing charge. The condition holds only within the light cylinder. The 
direction of particles in the lab frame is given by 



where 



/3 = /3[ | B + nxr/c, 



/?[, = -B ■ (O x r/c) + {[B • (n x r/c)] 2 - (fix r/c) 2 + 1} 1/2 , 



(3) 



(4) 



and B is the unit vector along the magnetic field in lab frame. The particle velocity is highly 
relativistic, so we have made the approximation |/3 | — > 1 in eq. (j3]). Thus the direction of radiation 
emitted tangential to the particle velocity vector is given by /3 in eq. ([3]). This direction is related 
to the periodic pulse phase. The observed phase eft is the sum of the azimuthal angle (p em at the 
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emission point and the relativistic time delay (jRomani fc Yadigaroglulll995l ): 



Rlc 



(5) 



where r em is the emission point, and Rlc the light cylinder radius. 



A certain mechanism is needed to fix the lower boundary of the particle acceleration region. 
In most works, including TCS08, the lower boundary is chosen as the surface of the last-open 
field lines of a rotating dipole in a vacuum. The field lines are calculated by eq. ([2]) from the 
neutron star surface, and the last-open ones are defined as being just tangential to the light cylinder 
and they form a magnetic surface from the polar cap. The numerical procedure is described by 
Cheng. Ruderman Zhanel (|2000l ) in detail and we follow it. In the outer gap model, if particle 
acceleration occurs in an open zone, the curvature radiation from the accelerated particles forms 
a narrow cone along the magnetic field lines in a frame rotating with angular velocity f2. These 
7-ray photons are converted by colliding X-ray photons to pairs, which tend to screen the 
accelerating electric field. However, there is no supply of pairs on the last-open field lines and 
hence no screening of the electric field, since the 7-ray photons are e mitted only toward higher 
altitudes above the last-open field lines (ICheng. Ho &; Rudermanlll986l). The 'real ' last-open field 



lines may be different from ones in a vacuum (jRomani 



1996 



Kojima &: Oogill2009l ). We therefore 



take into account this possible deviation of the boundary. We assume that dipole magnetic field is 
an approximation within the light cylinder and use eq. ([2]) as the global magnetic field structure. 
Even if the overall structure is not different so much, critical value between open and closed field 
lines is very sensitive to the boundary value at the surface. The "real" last-open field lines do 
not in general agree with those in vacuum field. Thus we introduce a parameter, altitude of the 
emission region as a correction factor in order to take into account the deviation of boundary from 
the vacuum field. In our model this parameter specifies the range of the emission region which 
is located above or below the last-open field lines within the light cylinder radius. Each different 
field line originating from the magnet ic polar region is parameteri z ed by magnetic colatitudes 9 m 
and azimuthal angles <f) m . Following ICheng. Ruderman Zhangl (|2000l ). we define open volume 
coordinates on the polar cap, (r ov , <j) m ), where r ov = m /6m'°{4>m)- The function is the 

magnetic colatitude of the conventional polar cap angle and generally depends on the magnetic 
azimuth cp m . The parameter r ov corresponds to the altitude of the emission region: The last-open 
field lines of a rotating di pole in a vacuum c orresp ond to r ov = 1, whereas those for higher altitudes 
have r ov < 1. Following iTakata fc Changl (2009j), the maximum value is chosen as r ov = 1.36 1//2 , 
which corresponds to the polar cap a ngle Om ~ 1.36 l / 2 0m' , obtained in the force-free limit by 



Contopoulos. Kazanas &: Fendtl (jl999l ) . We found that no significant caustics are formed in the sky 
map, even if the maximum value of r ov is increased. 

We assume that the radiation of different energy bands is emitted from different field lines 
characterized by altitude. The field line relevant to the 7- and X-rays is approximated as being 
the same one. The direction of the emission is tangential to the lines, and inward and outward 
directions are possible. Both location and direction affect the light curve profile of the energy 
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bands. Following the model by TCS08, the 7-ray radiation above 100 MeV is emitted by particles 
moving in an outward direction, whereas radiation at lower energy bands is emitted by those moving 
in both outward and inward directions. We use two conditions to constrain the emission region. 
First condition is the radial extension of the emission region. The outward emission is restricted 
to radial distances r n < Rlc-, and the inward one is restricted to r s < r < min(3r n , Rlc)- The 
outer boundary 3r n for inward emission comes from the results of dynamic model (TCS08), in 
which very few ingoing pairs are produced beyond the radial distance r > 3r n . Second condition 
is the azimuthal extension of the emission region. We use the magnetic azimuthal angle of the 
footprint of field line (i.e., the point where magnetic field line penetrates the neutron star surface) 
to characterize the field line for given v ov . Radial distance to the null charge surface on the field 
lines significantly depends on the magnetic azimuthal angle. In the outer gap model, most of the 
pairs are created around the null surface (TCS08). We expect that the gap activity is related to 
the distance to the null surface. Although the current density should be determined by global 
conditions, there is no study of the three-dimensional magnetosphere of an inclined rotator. In this 
paper, we assume that the field lines of both outward and inward emission are active only if the 
radial distance to null surface r n is shorter than Rlc- The azimuthal constraint is automatically 
satisfied for outward emission because the radial extension gives r n < Rlc- However, for inward 
emission the condition becomes strong. The radial extension r s < r < min(3r„, Rlc) allows for the 
regions r < Rlc on the field lines with r n > Rlc- They are not active, so that the corresponding 
regions should be excluded. The critical value 3r„ was obtained by fitting to Vela pulsar (TCS08). 
It is not straightforward to apply it to other sources. The mean free path A(r) of the pair creation 
process between the 7-ray and th ermal X-ray emi ssions from the stellar surface is estimated as 



A(r) ~ 5.6P 13 / 21 (-B s /10 12 G)" 2 / 7 r (|Tang et all 12008 ). The value at the null pointis A(r n ) found to 



be in a range of (2-3)r n for our samples. Our light curves especially peak positions are not changed 
even by adopting 2r n as the outer boundary for inward emission. 

Spatial distribution of the emissivity is approximated by the step function-type, but the peak 
positions weakly depend on the detailed emissivity distribution. 

We assume that the overall structure of the light curve comes not from the emissivity distri- 
bution, but from a bunch of many field lines in the observation, that is, caustics. The appearance 
of caustics strongly depends on the observational viewing angle £ and the intensity distribution. In 
this paper, we focus on the peak phases of the light curve, so we adopt a simple, uniform emissivity 
along all magnetic field lines, which is independent of both the magnetic azimuthal angle <ft m and 
the altitude r ov . The fitting does not completely reproduce the observations so, in Section 4, a 
simple improvement to the emissivity distribution is considered which leads to a much better fit. 

We now explain our fitting method. For fixed inclination angle a and viewing angle £, the 
light curve as a function of phase <j) depends only on the altitude r ov . The intensity is calculated 
in the range r ov < 1.36 1/2 with a bin width of 0.02. There are no significant caustics for large t ov . 

In the observed light curve, the reference phase (p = is assumed to be located at the radio 
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emission peak maximum in most studies (e.g., lAbdo et al.l (|2010al )). However in the model light 
curve, the conventional reference phase eft = occurs when the magnetic axis, spin axis and Earth 
line of sight lie all in the same plane. These two reference phases do not agree with each other 
since it is generall y assumed that radio em issions arise at non-zero altitude in most empirical 
studies. Following iRomani fc Wattersl (|2010l ). we allow a shift by —0.1 < 5(j) < 0.1 in the model 
reference rotation phase. This degree of freedom does not significantly affect the determination of 
the altitude parameter r ov , because we use the peak separation in the observed 7-ray and X-ray 
light curves which are emitted at the same v ov . For the sources showing a double-peak pulse shape 
in the observed 7-ray light curve, we use the peak separation. For those showing a single-peak, we 
use the phase separation between the 7-ray peak and one of X-ray peaks. This is the benefit of 
considering 7-ray and X-ray light curves simultaneously. Subsequently, we look for the altitude of 
the UV/optical emission region using the 7-ray upper limit of r ov . 



3. RESULTS 



In this section, we compare our model with pulse profiles observed at multiple wavelengths 
for seven pulsars. The sources are chosen using two criteria. One is that non-thermal pulses are 
detected in addition to the 7-ray and radio bands. Our concern is to explore whether or not the 
emission region for different energy bands is explained by the TCS08 model. The second criterion 
is that the geometrical parameters, a and £ are observationally cons trained by the relativistic 
Dop pler-boosted X-ray pulsar w ind nebula (PWN) (|Ng Romanill2008l ) or radio polarization data 
(e.g. iLvne Manchester! (|1988l )). The torus fitting method constrains the viewing angle £ only. A 
small allowed range of \a — £| < 10° is assumed for samples in which only £ is constrained due to 
the fact that radio emission from the pulsar polar region is detected. The geometrical parameters 
for the pulsars are listed in Table 1. We use these values, although there are some uncertainties in 
them. The results are summarized in Figs. [T] and [2j Fig. Q] shows the intensity map for outward 
(upper panel) and inward (lower panel) emission as a function of the altitude of the emission region 
r ov and rotational phase. The upper and lower panels in Fig. [2] are their pulse profiles for outward 
and inward emission in 7-ray and X-ray emission regions. 



3.1. Vela pulsar (PSR J0835-4510) 



We start with the Vela pulsar, which has been well studied to test the validity of our simple 
model. TCS08 considered this source, but they used geometric parameters that are slightly different 
from ours. 

Pulse profil es have been detect ed in optical to 7-ray bands. The observed pulse profile in 7-ray 
band by Fermi (jAbdo et al.ll2010bl ) shows a prominent double-peak structure and bridge emission 
between the two peaks. The first and second peaks are located at the phases <fi ~ 0.13 and <fi ~ 0.56, 
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respectively, and the separation is A(ft = 0.43. We show the intensity map for outward emission as 
a function of the altitude of the emission region and rotation phase in upper panel of Fig. [2(A). 
The emission altitude producing a peak separation A(ft = 0.43 is r ov ~ 1.05-1.06. 



The X-ray data from RXTE ([Harding et al.ll2002l ) also shows a double-peak structure but the 
second peak broadens toward early phase. The calculated intensity map is shown for outward and 
inward emissions in the upper and lower panels of Fig. [2(A). The main double peaks are located 
at the same phases as those in the 7-ray band, so they are interpreted as being formed by outward 
emission. The broad component before the second peak at (ft ~ 0.47 is associated with the caustic 
formed by the inward emission, as shown in the lower panel. We attempted a fit without the inward 
emission, but found that the inward emission is needed in order to reproduce the observed X-ray 
pulse profile. The necessity of inward emission was discussed in TCS08. Thus, the peak positions 
of 7- and X-ray pulses can be explained with the same value of r ol 



map, however, shows a minor peak at 
observed in the 7-ray band. 



■ 1.05-1.06. The contour 
0.8 formed by the outward emission. The peak was not 



We compare our m odel with the UV data of lRomani. Kargaltsev &: Pavlovl ([20051 ) and optical 
data of iGouiffes! (|l998l ). The pulse profiles in both bands are very similar, that is, they have a 
double-peak structure at the same phases. The peak phases however differ from those of the 7- 
and X-ray bands. The first peak of the UV/optical bands shifts to a later phase (ft ~ 0.27 and the 
second peak shifts to an earlier phase (ft ~ 0.46, so that the peak separation becomes smaller. It 
can be seen from Fig. [2(A) that such a double-peak structure corresponds to r ov ~ 0.65-0.80 for the 
outward emission. The corresponding inward emission cannot be detected since its observable range 
is v ov ^ 0.9, as shown in the lower panel. The choice of outward emission with v ov ~ 0.65-0.80 is also 
supported by the fact that the second peak at (ft ~ 0.46 in the UV/optical ranges is sharper than 
the first one at (ft ~ 0.27, because of their different dependence on t ov . Thus, we have reproduced 
the pulse profiles of optical to 7-ray bands by the caustics model without any detailed assumptions 
about emissivity. From the fitting model, we found three conditions for the emission region. (1) 
The UV/optical emission region is located at an altitude above the 7- and X-ray emission region 
of r ov ~ 1.05-1.06. (2) There is a separation of altitude between the X-ray and optically dominant 
emission regions. (3) The UV/optical emission range, Ar ov ~ 0.15, is broader than that of 7/X-ray 
emission regions, Ar ov ~ 0.02. 



3.2. PSR J0659+1414 

The pulsar PSR J06 59+1414 has also been observed in the 7-ray to optical ban ds. We use the 



7-ray data from Ferm i dWeltevrede et al J 120101 ). X-ray data from XMM-Newto n (iDe Luca et al 



20051 ). UV data from IShibanov et all (120051 ) and optical data from iKern et al.l (120031 ). The X- 
ray data is a combination of thermal (blackbody) and n on-thermal (power-law) em issions and is 



consistent with a cooling middle-aged neutron star (e.g., iBecker k, Trump erl (jl997l )). At soft X- 



rays, the pulse fraction is low and the pulsations are sinusoidal, as is typical for thermal emissions 
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from the surface of a neutron star with non-uniform temperature distribution. At higher energies 
(>1.5keV), where the non-thermal component dominates, the pulsed fraction increases and the 
profile becomes single peaked. We, therefore, consider the pulse profiles of hard X-rays (>1.5keV) 
only. 

The pulse profile in the 7-ray band shows a relatively broad single peak, which lags the radio 
maximum peak by 4> ~0.2 in phase. The peak in the non-thermal X-ray pulse is at <f> ~ 0.7-0.8, 
which is different from the phase of the 7-ray peak. This phase difference cannot be ignored, 
although the peaks in the 7- and X-ray data are rather broad, and hence the difference may be 
diminished somewhat by including the phase error. We interpret these pulse profiles as being 
emissions at different phases and different directions: the peak of 7-ray is formed by outward 
emission, whereas that of X-ray is formed by inward emission. The intensity maps are given in 
the upper panel (outward emission) and lower panel (inward emission) of Fig. QjB), respectively. 
From this figure, we see that a peak separation Acp = 0.55 between the 7-ray and X-ray data 
corresponds to an emission altitude of r ov ~ 1.13-1.14 by shifting the reference phase by 5(f) = 0.06 
to an earlier phase. The emission altitude cannot be fixed without the X-ray data: a shift of peak 
phase is allowed, so r ov is unknown. This ambiguity is removed by considering multi-wavelength 
light curves. From the intensity map, we expect another very shar p caustic at <t> ~ 0.65, but there 



is no counterpart in the 7-ray observations (top panel of Fig. 4 in IWeltevrede et al.l (2010)) . We 



discus s this missing peak in Section 4.2. In the lower right panel of Fig. 4 in IWeltevrede et al 
(|2010l ). the light curve for inward emission is given. The X-ray data may be a comb ination of 
inward and outward emissions, but the X-ray profile observed by iDe Luca et al.l (J2005J) is similar 
to that of inward emission only. This means phenomenologically that outward emission of X-rays 
is weak in this source. 



The pulse profiles in the UV (jShibanov et alj|2005l ) and optical bands (jKern et alj 12003 ) are 
almost the same shape and have a clear double-peak structure unlike the single peak in the 7- 
and X-ray bands. The first peak at eft ~ 0.02 is later than the 7-ray peak and the second peak at 
(f> ~ 0.10 is later phase than the X-ray peak. From the upper panel of Fig. QjB), we see that the 
observed first peak can be reproduced by outward emission at an altitude r ov < 1.10. Here we have 
a weak condition because the peak position depends only weakly on t ov . For the second peak, the 
inward emission forms caustics for 0.90 < r ov < 1.04. Thus, the altitude of the emission region in 
the UV/optical bands is identified as v ov ~ 0.90-1.04, where the lower limit is set by a coarse bin 
of the phase in the observational data. 



Kern et al.l (|2003 ) have already investigated the multi-wavelength light curve of this pulsar 



using a similar method to ours, but could not e xplain the profile using ge ometrical parameters 
which are consistent with radio polarization data ([Everett fc Weisbergil200ll ). The reason for this 
is that the lower boundary of the emission region was chosen as the last-open field lines in the 
vacuum dipole field, that is, r ov = 1.0. In our analysis, by allowing r ov > 1.0, the phase of peaks 
can be successfully fitted by using observed geometrical parameters. This suggests that the actual 
lower boundary of the gap is slightly different from the last-open field lines in a vacuum dipole. 
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3.3. PSR J0205+6449 



The 7-ray pulse profile observed by Fermi (jAbdo et al.ll2009bl ) shows a double-peak structure. 
The first peak is offset from the radio peak by (ft ~ 0.08, and the second is at (ft ~ 0.57. The 
separation is A(ft = 0.49. The X- ray data from RXT E ([Livingstone et al.ll2009i ) are consistent with 
the results from XMM-Newton (IKuiper et al.ll2010l ). The spectral shape can be fitted by a power 
law such that m ost of the emission is non-thermal and the thermal component is constrained by 
the upper limit (jKuiper et al.ll2010l ). The observed X-ray pulse profile shows two peaks aligned in 
phase over a wide energy range of ~0. 5-270 keV, and is also very similar to that of the 7-ray band. 

We show the intensity maps for outward and inward emission in the upper and lower panels of 
Fig- QIC), respectively. As seen in the upper panel, the emission altitude for the double-peak with 
A(ft = 0.49 is r ov ~ 0.97-0.98. A shift of the reference phase 5 (ft is not necessary in this source. As 
argued in TCS08, outward emission dominates in the light curve for a young pulsar with a strong 
non-thermal X-ray component, like the Crab pulsar. PSR J0205+6449 is the youngest pulsar in 
our sample (characteristic age r c ~ 5 x 10 3 yr) and shows rather strong non-thermal radiation in 
the X-ray band. Thus, it is likely that only outward omission contributes to the observed X-ray 
light curve in this pulsar. 



3.4. PSR J2229+6114 



The light curve for 7-rays observed by Fermi (jAbdo et al.ll2009d ) shows an asymmetric single 
peak at (ft ~ 0.49. The tail extends down to (ft ~ 0.2. The peak position depends slightly on the 
energy range above 100 M eV, but the amoun t of shift is only ~ 0.04. The X-ray pulse profile 



observed by XMM-Newton (jAbdo et al.ll2009d ) shows a double-peak structure. No peak is seen at 
the 7-ray peak phase. The separation between first X-ray peak and the 7-ray peak is A(ft ~ 0.32, 
and the separation between the second X-ray peak and the 7-ray peak is A(ft ~ 0.14. 

The intensity map for outward emission is shown in the upper panel of Fig. QJD). We consider 
the formation of the 7-ray peak and two X-ray peaks as being due to outward emission only. Such 
a solution is possible by choosing an emission altitude of r ov ~ 1.01-1.02 with a small phase shift 
5(ft = 0.03. The intensity map for inward emission shows a sudden decrease in the number counts 
for r ov < 1.06. The peak becomes broad and hence the contribution of the inward emission is not 
very important. We show the light curve of inward and outward emission for r ov ~ 1.01-1.02 in 
Fig. DJD). The outward emission curve is very similar to the observations in the 7- and X-ray 
bands. 

In this pulsar, we need to use all the light curves simultaneously in order to determine the 
range of r ov . Since emissions with smaller values of r ov are not seen, if 7-ray and optical emission 
regions are separated Ar ov > 0.10, similar to the Vela pulsar and PSR J0659+1414, we predict that 
an optical pulse profile cannot be observed or will be only very weakly detected. This is consistent 
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with the fact there have been no reports of the detection of a pulse in the lower energy band for 
this pulsar. 



3.5. PSR J1420-6048 



The 7-ray light curve from Fermi (jWeltevrede et al.ll2010l ) shows a broad peak at (ft ~ 0.2-0.5. 
This peak may consist of two components, but it is not clear in the current photon statistics. 
The X-ray pulse profile from ASCA is detected weakly at a marginal level, and shows a broad 
peak at <ft ~ 0.6-0.7, which is different from the 7-ray peak (Roberts. Romani &: Johnston! l200ll ; 



Becker 2009 ). Recently, in the table of Marelli. De Luca Sz Caraveol ( 2011 ) they list this object 



as 

non-thermal dominated source in X-ray. The pulsed X-ray profile is likely to originate from the 
non-thermal component. 

Since the light curve of this pulsar and its geometrical parameters are similar to those of 
PSR J0659+1414, we adopt the same interpretation. That is, the 7-ray peak is formed by outward 
emission, whereas the X-ray peak is formed by inward emission. From the intensity map in Fig.QjE), 
an emission altitude of r ov ~ 1.10-1.11 corresponds to one broad peak at (ft ~ 0.2 — 0.5 by outward 
emission and another at (ft ~ 0.6 — 0.7 by inward emission. Here a small shift 5(f) = 0.10 toward 
earlier phase is used. Since there are similarities in both the 7- and X-ray light curves and the 
geometrical parameters between this pulsar and PSR J0659+1414, we expect a similar double-peak 
pulse profile in the optical band, if it is detected. 



3.6. PSR J2021+3651 



Observations in the 7-ray band have been obtained by Fermi (jAbdo et al.ll2009dl ) and AGILE 
(|Halpern et al.l 12008). The observed light curve shows a sharp double-peak structure. The first 
peak is offset from t he radio peak by (ft ~ 0.16 and the two peaks are separated by Acft ~ 0.47. The 



X-ray light curve in lAbdo et al.l (|2009dl ) shows a relatively sharp peak associated with first peak 
in the 7-ray light curve albeit with weak photon statistics. In this paper, we assume that at least 
the first peak is n on-thermal in origin . Th e possible contribution of non-the rmal X-ray emissions 
is also discussed in lHessels et al.l (|2004l ) and lVan Etten. Romani Na (|2008l ) . We expect that this 
assumption will be tested by phase-resolved spectra from future observations. 

As seen in the upper panel of Fig. HHF), the emission altitude 0.97-0.98, for which there 

is a 7-ray double-peak profile with separation = 0.47 and a relative shift 5(p = 0.06 toward later 
phase. The peak of non-thermal emission at (ft = 0.15-0.20 in the X-ray light curve is found to be 
formed by outward emission only. The relatively weak second peak in the X-ray band is consistent 
with the case of PSR J0205+6449. Thus, the three model parameters for this pulsar and PSR 
J0205+6449 are very similar, as shown in Table 1. 
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3.7. PSR J1057-5226 



The 7-ray light curve from Fermi (jAbdo et al.ll2010d ) s hows a broad p e ak at (p ~ 0.25-0.65. 
This probably consists of two components, but it is not clear. |Pe Luca et al.l (|2005l ) extract only a 
power-law component of the X-ray light curve and their Fig. 13 shows a two peaks at <f> ~ 0.2-0.3 
and 4> ~ 0.9-1.0, although the data are very coarse. We regard the light curve as being produced 
by a non-thermal X-ray component. 

The observed light curve may be regarded either as a broad peak consisting of weak peaks 
and a relatively bright bridge emission or as a result of the range of the emission region widening 
towards lower altitudes. In the latter interpretation the fitted r ov is only a lower limit. We thus 
focus on the width of the 7-ray peak. From the upper panel of Fig. QJG), the emission altitude 
is r ov ~ 0.93-0.94 Even if we assume a double-peak struct ure with the first peak at (j> ~ 0.31 and 



the second peak at (j) ~ 0.59 following lAbdo et al.l (|2010d ). we have r ov ~0. 90-0. 91, which is very 
similar to the value obtained above. Thus we have r ov ~0. 90-0. 95 in either case. The phase shift is 
5<f) = 0.10 in this pulsar. The first peak in the X-ray light curve is formed by outward emission, but 
the second one cannot be produced for the same altitude. This may be a drawback to our model, 
but the present X-ray data are coarse and a much more precise non-thermal X-ray light curve is 
needed. 



4. DISCUSSION 

4.1. Statistical properties of the emission region 

In the previous section, we have shown that the peak phases of seven pulsars emitting 7- and 
X-rays can be successfully fitted using the TCS08 outer gap model, in which both 7-rays and X-rays 
originate from the same magnetic field line characterized by an altitude r ov . The parameter r ov > 1 
is needed in the light curve fitting for some sources. Moreover, the inclusion of inward emission for 
X-rays causes a variety of pulse profiles in both bands. The parameter r ov could not be determined 
solely using 7-ray data for a single 7-ray peak pulsar. But, by considering the X-ray light curve, 
the parameter is uniquely determined for PSRs J0659+1414, J2229+6114 and J1420-6048. 

It is worthwhile to explore the general dependence of the altitude r ov on other characteristics 
if any, although there may not be enough data for a proper statistical analysis. In Fig. [3l r ov is 
plotted as a function of inclination angle a, spin-down luminosity L s a, characteristic age r c and 
surface dipole magnetic field B s . We found that there is a significant correlation between r ov and 
the inclination angle a only; the relations of r ov with the other parameters are very weak. This 
correlation suggests that the deviation from a vacuum rotating dipole field is large for small incli- 
nation angle. It is yery in teresting to compare this result with that in a force-free magnetosphere. 



Bai Spitkovskvl (|2010bl ) proposed that the separatrix layer at an altitude of 0.90-0.95 times the 



height of the last-open field line is relevant to emissions in a three-dimensional inclined force-free 
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magnetosphere. This altitude, which is not exactly symmetric with respect to the magnetic az- 
imuthal angle (p m , but can be approximated by the value at (fi m = 0, is plotted in Fig. [3] as purple 
downward and blue upward triangles. Two linear fitting lines are also shown. The altitude r ov de- 
creases with the inclination angle a in both our model and the separatrix layer model of a force- free 
magnetosphere. However, the emission region in the separatrix layer model extends even outside 
the light-cylinder, whereas ours is well localized around null points. Accounting for this difference 
may be important for further improvement of the model of the emission region based on a force- free 
magnetosphere. 

The thickness of the gap r egion, w, is not known, but it is sometimes ass umed to decrease with 



the spin-down luminosity L$d (jWatters et al.ll2009l ; lRomani Watters!l2010l ). We have w = 1— r, 



if the lower boundary of the gap is fixed as the last-open field line in the vacuum dipole field. 
This assumption is tested in the lower left panel of Fig. [3l in which the relation (1 — r ov ) ~ 
(Lso/W^ergs- 1 )- 1 / 2 is plotted as a light green curve. (The curve is not fitted to the data 
points.) This suggests that the assumption of maximum altitude, r ov = 1.0, is not a good one. 
This discovery affects expected number of the 7-ray pulsars in the observation. From geometrical 
reason, the pulsed emission by caustics is limited to a certain range between inclination and viewing 
angles. 



Romani &; Wattersl (|2010l ) showed the range of observable pulsars with r ov =0.95, 0.90 and 
0.70 for outer gap model in their Fig. 16. We recalculate it and show the result in Fig. HI The 
observable range of viewing angle £ is below the curves. Our finding in Fig. [3]is that r ov is a function 
of the inclination angle, which is similar to that of the separatrix layer model. We also show the 
observable range by the empirical relation obtained in Fig. [3] as black solid line, for which the 
altitude is chosen as 0.925 times the height of the last-open field line in force-free magnetosphere. 
The figure shows that sources with low inclination and viewing angles become observable. For 
example, pulsar with the inclination angle a = 30° can be detected for £ > 60° for r ov =0.95, 
but for £ > 30°. Thus expected number increases approximately twice for sources with the low 
inclination and viewing angles. 



4.2. The phenomenological limitation for emissivity 



The caustic model considered in Section 3 provides peak positions consistent with observation, 
but there are also some additional, unseen peaks. These are interpreted as being prohibited by 
some mechanism. In this section, we consider an improvement to our model that takes into account 
a very simple distribution for the emissivity. Detectable 7-rays are radiated with large multiplicity 
by the pair plasma in the gap region. Therefore, the mean free path of a 7-ray photon should be 
less than light cylinder radius ([Takata &: Chans 120071) . The pair creation mean free path is given 
by A(r) ~ 5.6P 13 / 21 (B s /10 12 G)~ 2 / 7 r (|Tang et al.l 120081 ) for an assumed limiting distance to the 
null point r nt n m . The position of the null point of inclined pulsars, where the accelerating electric 
field arises, depends significantly on magnetic azimuthal angle, so the intensity of 7-ray emission 
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also depends on the magnetic azimuthal angle. Therefore, active field lines should be limited in the 
azimuthal direction. By taking into account the azimuthal extensions with A(r n ^ TO ) 0.2 — 0.7i?^c 
listed in Table 1, the fits of the resultant light curves, which are shown in Fig. [5j become better. 
In the same figure, we also show the radial distance to the emission points of the observed photons 
against the rotation phase. Note that, for the Vela pulsar, the minor third peak at <p ~ 0.8 still 
remains even after the inclusion of the azimuthal extension limit. The corresponding radial distance 
of emission points is relatively large, so that the photon energy is expected to be soft. The third 
peak is not observed in the GeV band, but may appear in a much lower energy band. At least, the 
minor third peak of the X-ray light curve appears to be associated with the same caustic. 

We have also tried to improve the X-ray light curve with some other simple assumptions, but 
have not had good results. The reason for this is that there are many ways for X-ray emitting par- 
ticles to be created: via thermal, magnetospheric emissions and magnetic pair creation. Therefore, 
the three-dimensional effect of the propagation of 7-ray photons and soft X-ray photons is very 
important. Without it we cannot successfully explain the light curve. 



4.3. The location of the UV/optical emission region 

We explored the UV/optical emission region for Vela and PSR J0659+1414. The results are 
qualitatively similar: the altitude range of UV/optical emission, Ar ov ~ 0.15, is broader than that 
for 7- and X-rays, Ar ov ~ 0.02; and both emission regions are not continuous and connected, but 
are widely separated. The separation may come from two competing mechanisms: a decrease of 
emissivity and an increase of synchrotron intensity in the UV/optical bands with altitude. 

As discussed in TCS08, outward emission is generally dominant in UV/optical emission, as 
shown in their Fig. 4. The explanation is the following. The number of created pairs is the main 
cause of the difference between inward and outward emissions in the UV/optical bands, since the 
collision angles with magnetospheric X-rays are not different for outgoing and ingoing 7-ray photons 
in the acceleration region. More outgoing 7-rays are emitted, and hence more outgoing secondary 
pairs are produced. Thus, the synchrotron emission for outgoing secondary pairs produced by 
magnetic X-rays is brighter than that for the ingoing secondary pairs. The observed flux strongly 
depends on the geometrical configuration. The outward emission in the UV/optical bands may not 
point toward us even though the intrinsic emission is strong. Our results show that the peaks in the 
Vela pulsar can be explained by outward emission alone, while those in PSR J0659+1414 require 
both inward and outward emissions. Our result suggests that outward emission is significantly 
suppressed in PSR J0659+1414, to the level of the intrinsically weak inward emissions. The stronger 
component is hidden, because the observable altitude range is narrow, as shown in Fig. QJB). This 
may explain the fact that UV/optical flux is sm aller than the value extrapolated from non-thermal 



X-rays, as seen Fig. 4 of iMignani et al.l (|2010l ). whereas the flux coincides with the extrapolation 
in the Vela pulsar. This interpretation may be tested in PSR J1057-5226. We also suggest that 
this difference is the reason why PSR J0659+1414, which has similar geometrical parameters to 



-14- 



the Vela pulsar, has an observable optical spectrum (jMignani. Pavlov fc Kargaltsevll201Cll ). and the 
flux is slightly smaller than the extrapolation from non-thermal X-ray emission. The pulse profile 
has not yet been determined, but the peaks should appear at a phase 0.3< <p < 0.6 and be due to 
outward emission. 



5. SUMMARY 

We have calculated the light curves of emissions using the TCS08 outer gap model and com- 
pared them with observed multi- wavelength light curves. We find that the model can successfully 
explain the peak positions of multi-wavelength light curves. In order to determine the altitude of 
the emission region, the observed X-ray light curve is important, especially when there is a single 
peak in the 7-ray light curve. 

The fit of a light curve based on a simple emissivity distribution can be improved by taking 
into account the limitation of azimuthal extension in which a reasonable value of the 7-ray photon 
mean free path is adopted. The resulting difference between model and observed 7-ray light curves 
becomes small; however, there may still be an unseen peak, such as the minor third peak in Vela. 

The best-fit values of the altitude of the emission region for PSRs J0659+1414 and J1420-6048, 
suggest a deviation from the last-open field lines of a vacuum dipole field. The real last-open field 
lines lie inside those of vacuum dipole field, r ov < 1.0. This shift suggests that the lower boundary is 
very similar to that of a force-free magnetosphere. We find that the altitude of the emission region 
is correlated with inclination angle. This relationship is also very similar to that in a force-free 
magnetosphere. The lower boundary of emission region has been assumed to r ov = 1 so far, but 
our model fits do not support it. This modification of the boundary of the magnetosphere suggests 
that the pulsars with low inclination and viewing angles are lik ely to be detectable. Thus th e 



expected number in the f uture observation in the previous works (jTakata. Wang &: Chena (|201ll ) 



Watters fc Romanil (|201lh ) is underestimated for the sources with low inclination and viewing 



angles. 
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TABLE 1 Pulsar parameters 



Name 


log(L SD ) 




log(B.) 


a 




Reference 


r OT (7-, X-ray) 


r M (UV/optical) 








(erg s _1 ) 


(kyr) 


(G) 


(degrees) 


(degrees) 








{Rlc) 


(Rlc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


J0835-4510 


36.84 


11 


12.53 


72 


64 


1,2 


1.05-1.06 


0.65-0.80 


0.25 


0.23 


J0659+1414 


34.58 


110 


12.67 


29 


38 


3 


1.13-1.14 


0.90-1.04 


0.30 


0.60 


J0205+6449 


37.43 


5 


12.56 


78 


88 


2 


0.97-0.98 




1.00 


0.71 


J2229+6114 


37.35 


11 


12.31 


55 


46 


2 


1.01-1.02 




0.40 


0.29 


J1420-6048 


37.00 


13 


12.38 


30 


35 


6,7 


1.10-1.11 




0.50 


0.42 


J2021+3651 


36.53 


17 


12.50 


75 


85 


5 


0.97-0.98 




0.40 


0.40 


J1057-5226 


34.48 


540 


12.03 


75 


69 


4 


0.93-0.94 




0.10 


0.20 



NOTES. -Col. (l):Pulsar name. Col.(2),(3) , (4):The spindown luminosity, the characteristic age and the strength of surface mag- 
netic field, which we adopt in lAbdo et al.1 ( 2010a ). Col.(5):The inclination angle. Col.(6):The viewing angle. Col. (7): Reference 
for cols. (5) and (6). Col.(8):The emission altitude in the 7-ray band. Col.(9):The emission altitude in the optical/UV 
band. Col.(10):Assumed limit of radial distance to null point. Col.(ll):Mean free path for 7- ray photons at th e limit 
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Fig. 1. — The intensity maps for seven pulsars. Upper left to right: PSRs J0835-4510 (A), 
J0659+1414 (B) and J0205+6449 (C). Lower left to right: PSRs J2229+6114 (D), J1420-6048 
(E), J2021+3651 (F) and J1057-5226 (G). In each sample, upper panel is outward emission and 
lower is inward emission. The blue horizontal arrows show best fit values of r ov for 7-ray and X-ray 
emission regions. The red vertical arrows show the phase of peaks. The red horizontal arrows in 
(E) and (G) show the phase range of broad peaks. 



-20- 



(A) (B) (C[ 

- — IT 1 1 n 1 3 ( 1 rn 1 





\ 

XJ 


i i 

VI 


1 




1 1 


1 


1 1 1 


H 




H h- 

u i i 1 


1 1 
1 1 


1 1 1 


1 


H — hr 


J: 

i "i 





i 1 




1 


I 1 

: pj 


i i 

1, , 




i i i 


1 

u 


1 


i i 


1 

l I 1 




u 

1 


1 


,1 , 1 


i i i i 




f 




-\ — I — h 

> 

*T 1 1 


1 


r 


i 


nl 1 1 1- 

1 f 
III 



0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 
<£ €> <i> <J> 

Fig. 2. — The calculated light curves for 7- and X-ray emission region. Upper left to right: PSRs 
J0835-4510 (A), J0659+1414 (B) and J0205+6449 (C). Lower left to right: PSRs J2229+6114 
(D), J1420-6048 (E), J2021+3651 (F) and J1057-5226 (G). In each sample, upper panel is outward 
emission and lower is inward emission. The vertical axis is in arbitrary units. The red and blue 
short-dashed vertical lines show the phase of 7-ray and X-ray peaks as in Fig [TJ The peaks of 
PSRs J1420-6048 and J1057-5226 are so broad that the phase range is within two red vertical lines. 
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Fig. 3. — The relation between r ov and inclination angle (upper), spin-down luminosity (lower 
left), characteristic age (lower middle) and surface magnetic field (lower right). The altitudes cor- 
responding to the separatrix layer model are shown as purple downward and blue upward triangles 
in the upper panel. The two lines are linear fitting lines for the separatrix layer model. The light 
green curve in the lower left panel shows the relation (1 — r ov ) = (10 33 ergs~ l / Lsd) 1 ^ 2 ■ 
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Fig. 4. — The observable range for 7-ray pulsars in the a-£ plane for the outer gap model. Black 
solid curve shows the boundary of observable pulsars using linear fitting line for the separatrix layer 
model. Red dashed, blue dash-dotted and purple dotted curves show the boundary with r ov =0.95, 
0.90 and 0.70, respectively. Light-green curcles show the pulsars in Table 1. 
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Fig. 5. — The distribution of the radial distance to the null point of the field line on which observed 
photons are emitted (each upper panel) and the light curves that are restricted by the azimuthal 
extension limit as a function of the rotation phase (each lower panel) for seven pulsars. The color 
shows the radial distance to the emitting point as 0.0 < r/R^c < 0.2 (red), 0.2 < r/Ric < 0.4 
(light green), 0.4 < r/R LC < 0.6 (blue), 0.6 < r/R LC < 0.8 (purple), 0.8 < r/R LC < 1.0 (light 
blue). The values r n n m for each pulsar are shown as black long-dashed horizontal lines. The 
red short-dashed vertical lines show the phases of the 7-ray peaks. For PSRs J1420-6048 and 
J1057-5226, the vertical lines show the phase range of broad peaks. 



